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Abstract 

We use a two-component model to describe the production mechanism of D mesons 
in vrp colHsions. The model combines the usual QCD processes plus fragmentation and 
recombination of charm that has been produced by nonperturbative QCD mechanisms. 
A hard charm component in the pion must be responsible of the particle anti-particle 
production asymmetries observed. 
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1 Introduction 



Recent measurements of charm meson production in irp collisions |1|, ^ 0, ^, Q indicate that 
there are important nonperturbative QCD phenomena in the production process. 
The so called "leading effect" or particle anti-particle production asymmetry is much bigger 
than predicted by next-to- leading (NLO), QCD predictions. 

In this letter we study the xp distribution of and mesons in the framework of a 
two-component model. The production of D mesons in the model is assumed to take place 
via two different processes, namely QCD parton fusion with the subsequent fragmentation of 
quarks in the final state and conventional recombination of valence and sea quarks present 
in the initial state. 

The asymmetry obtained with the conventional soft charm component does not reproduce 
the experimental results. We show that with the presence of a hard charm component in 
the pion, the sum of these two processes gives rise to an enhancement at large values in the 
xp distribution of leading mesons. 

Unlike the |]^, |^ where the asymmetry is generated by the presence of the ud diquark in 
the initial hadron, the meson production asymmetry seems to imply the presence of a hard 
charm component in the pion. 

To quantify the difference in the production of leading and non-leading particles, an asym- 
metry A is defined 

^ a{leading) — a{nonleading) 
a{leading) + a{nonleading) 

In the reaction TT{du) - proton{uud) the D^{cu) and D~{cd) are leading mesons while 5° 
and are non-leading. 

Several models have been proposed although none fully account for the difference in pro- 
duction between leading and non-leading particles. The intrinsic charm model [§] postulates 
the existence of quantum fluctuations in the beam particle that bring about Fock states 
containing cc pairs. The cc quarks having the same velocity as the original valence quarks 
are likely to coalesce forming leading particles. Although the shape of the asymmetry as 
a function of xp predicted by this model is similar to the experimental measurement, the 
prediction is to low for the whole xp interval. 

In it is pointed out that the annihilation of the u from the proton and the u in the pion 
would liberate the d of the pion which can then recombine to form a D~ and certainly not 
a D^. The will not be enhanced because it cannot be formed with this simple annihi- 
lation diagram. The recent measurement of significant asymmetry in Ac production in pp 
collisions however, shows that there must be other production mechanisms at work that 
must account for the asymmetry. 

In two component models ||^, |T0[ the production of charm mesons by parton fusion is the 
same for and D^. In section 2 the formalism to obtain the cross section by parton fusion 
is shown for charged and neutral charm mesons. The recombination of charm as the second 
component is discussed in section 3. Recombination, is different for the charged mesons 
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where a parton process favors the formation of D enhancing the D asymmetry. This will 
be discussed in section 4. 



2 Parton fusion production of charmed mesons 

In the parton fusion mechanism the D^{D'^) mesons are produced via the qq{gg) — > cc 
with the subsequent fragmentation of the c(c) quark. The inclusive xp distribution of the 
charm mesons in irp collisions, is given by HIT 



daPf 



y/s J Hab{xa,Xb, Q'^)^ ^^^^ dzdp'^dy , 



(2) 



where 



HabiXai Xbi Q 



^a,b {qa{Xa,Q )qb{Xb,Q 



+ qa{Xa,Q'^)qb{Xb,Q'^)) i \qq 

' dt 

+ga{Xa,Q'^)gb{,Xb,Q'^)—2 \gg, 

dt 



(3) 



and Xa and Xb being the parton momentum fractions, q{x, Q"^) and g{x, Q"^) the quark and 
gluon distribution in colliding hadrons, E the energy of the produced c-quark and -Dd/c (z) 
the appropriated fragmentation function. In eq. |^, is the squared transverse momentum 
of the produced c-quark, y is the rapidity of the c quark and z = xp/xc is the momentum 
fraction of the charm quark carried by the D. The sum in eq. |]runs over a,b = u, u, d, d, s, s. 
We use the LO results for the elementary cross- sect ions ^ 



I,, and f 



199- 



da 

di 



ira: 



{Q"^) cosh (Ay) + ml/rh 



9mi 



[1 + cosh {Ay)Y 



(4) 



da 
di 



Tca: 



(Q2) 8cosh{Ay) 



\99~ 



96m4 



cosh (Ay) + + 



2mt 



nit 



(5) 



■ml + p\.. 



[1 + cosh {Ay)] 

where Ay is the rapidity gap between the produced c and c quarks and m\ 

In order to be consistent with the LO calculation of the elementary cross sections, we use 
the GRV-LO parton distribution functions [^, and apply a global factor K ~ 2 — 3 in eq. 
^ to take into account NLO contributions [IB . 

We take = 1.5 GeV for the c-quark mass and fix the scale of the interaction at 
Following |TD[, we use a delta fragmentation function Dpi/^^z 
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2ml 
which seems 



to describe experimental data better than the Peterson fragmentation function. 
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3 Charmed meson production by recombination 



Sometime ago V. Barger et al. explained the spectrum enhancement at high in Ac 
production assuming a hard momentum distribution of charm in the proton. According with 
them charm anti-charm pairs which give rise to the flavor excitation diagrams (see fig. 1) 
are not intrinsic but generated by QCD evolution of the structure functions. 
In this framework, following the gluon scattering process of fig. 1 the charm quarks will 
fragment into charm hadrons. When the c is scattered, the spectator quark could recombine 
with the u valence quark of the pion to form a D^{cu) or less frequently with an antiquark 
from the pion's sea. On the other hand, when the c quark is scattered, the spectator quark 
c could recombine with d valence quark from the pion to form a D~{cd) meson. 



As pointed out in the charm hadrons resulting from the scattered charm quark, popu- 
late the low xp region of the cross section, while those originating from the spectator quark 
dominate at high xp- 

Here we assume a QCD evolved charm distribution, of the form proposed by V. Barger 
et al. m 

xc{x, {Q^)) = Nx\l - x)\ (6) 

with a normalization N fixed to 

' dx ■ xc{x) = 0.005 (7) 



and I = k = 1. With this values for / and k one tries to resemble the distribution of valence 
quarks. In contrast with the parton fusion calculation, in which the scale of the inter- 
action is fixed at the vertices of the appropriated Feynman diagrams, in recombination the 
value of the parameter should be used to give adequately the content of the recombining 
quarks in the initial hadron. We used = Arn^ and therefore, the integrated charm quark 
distribution would take the value given in eq.7. 

The production of leading mesons at low pt was described by recombination of quarks 
long time ago . 

In recombination models it is assumed that the outgoing hadron is produced in the beam 
fragmentation region through the recombination of the maximum number of valence quarks 
and the minimum number of sea quarks of the incoming hadron. The invariant inclusive Xp 
distribution for leading mesons is given by 

2E da'''''' f^'F dxidx2 ^ , , „ , , , , 

—= =/ F2[Xi,X2)R2[Xi,X2,Xf) (8) 

a/SCT dXp Jo Xi X2 

where xi, X2 are the momentum fractions and -F2(xi,a;2,) is the two-quark distribution 
function of the incident hadron. R2{xi,X2,Xf) is the two-quark recombination function. 
The two-quark distribution function is parametrized in terms of the single quark distributions 

F2 (Xi, X2, ) = l3Fd^^al {Xl) Fc^sea {X2) (l - Xi - X2) , (9) 
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with Fq (xj) = Xiq (xj). We use the GRV-LO parametrization for the single quark distribu- 
tions in eqs. ^ The charm contribution however, is parametrized with the distribution of 
eq. 6. It must be noted that since the GRV-LO distributions are functions of x and Q^, then 
our F2 {xi^X2) also depends on Q^. The recombination function is given by 

R2 {xd, Xc) = a^^^^6 {xd + Xc - xp) , (10) 

X p 

with a fixed by the condition dxp{l/a)da'^^'^/dxp = 1. 



4 and D^{D^) total production 

The inclusive production cross section of the D is obtained by adding the contribution of 
recombination eq. |^ to the QCD processes of eq. ^ 

_ daPf ^ da"^^^ 
dxp dxp dxp 

The resulting inclusive D production cross section da^°^/dxp is used then to construct the 
asymmetry defined in eq. [^. 

In the vrp reaction however, the recombination of a D^{dc) gets a higher probability than 
the recombination of a uc 0. 

The process that releases a d quark from the incident vr to form a D~{dc) ( after recombina- 
tion with a c-quark from the sea ) is also present in the formation of a D^. In this case the 
released quark should be u. In the D~ case however, there is an additional mechanism that 
releases the d quark from the pion, namely the fusion of a w-quark from the proton with the 
u from the pion. This is correctly incorporated by the fact that a higher contribution from 
the recombination is needed to describe the D~ asymmetry while D^{D^) production and 
its production asymmetry requires a small contribution from the recombination process. 
Fig. 2 shows the model prediction for the D~ and cross section and experimental results 
from the WA82 Collaboration 0. The corresponding production asymmetry is shown in fig. 
3 together with all the experimental measurements available. Fig. 4 shows the cross section 
for and mesons produced in 7[~p collisions and the experimental points from the 
WA92 Experiment The corresponding production asymmetry for and Z)° is shown 
in fig. 5. By looking at the fit of the cross section in fig. 4 one can see that the asymmetry 
obtained for D^{D^) could actually change drastically by slightly changing the curves. More 
statistics is needed to improve the confidence of the predicted asymmetry. 

5 Conclusions 

In an earlier work the production asymmetry of Ac was described using the same recom- 
bination scheme used here. The presence of a diquark in the initial state, plays an important 
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role in Ac production. However, it is not possible to describe the asymmetry in charm mesons 
using the GRV structure functions for the charm in the pion. A hard charm component must 
be present. This fact has been pointed out before 0] but the asymmetry in the intrin- 
sic charm model M does not seem to fit experimental results very well. The hard charm 



component proposed in and a recombination scheme, give a good description of particle 



anti-particle production asymmetries for mesons. 
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Figure Captions 



Fig. 1: Heavy flavor excitation diagram. This and similar diagrams are included in a NLO 
QCD calculation. However, in oder to reproduce particle anti-particle asymetries a 
hard charm component in the pion is needed. 

Fig. 2: Production cross section for D~ and mesons in yr^p collisions. The dashed line 
represents the parton fusion contribution with a delta fragmentation function. The 
dotted line gives the contribution from recombination for D~ . The solid line is the 
sum of the two contributions. 

Fig. 3: Measured production asymmetry for D~ and and the model incorporating a hard 
charm component in the pion (solid line). The dotted line shows the intrinsic charm 
model prediction. 

Fig. 4: Production cross section for and mesons produced in vr^p collisions. Experi- 
mental points p and two component model prediction (solid line). 

Fig. 5: Measured production asymmetry for and and the model incorporating a hard 
charm component in the pion. The horizontal line at A{xp) = is for reference only, 
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Figure 3: 
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Figure 5: 
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